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Abstract 

■ The relativistic structure of the self-energy of a nucleon in nuclear matter 

£S) \ is investigated including the imaginary and real components which arise from 

the terms of first and second order in the NN interaction. A parameterized 
form of Brueckner G matrix is used for the NN interaction. The effects of the 

, terms beyond the DBHF approximation on quasiparticle energies and the optical 

potential for nucleon- nucleus scattering are discussed. 

Ij 1 Introduction 

During the last few years the attempts to derive the ground-state properties of nuclear 
systems from a realistic nucleon-nucleon (NN) interaction, have been promoted very 
much by the understanding that relativistic effects may be non-negligible in such in- 
vestigations. These ideas were originally developed within the various versions of the 
phenomenological Walecka model|l[]. The Dirac structure of the NN interaction with a 
strong repulsive component originating from the exchange of the u vector meson and an 
attractive component of medium range described in terms of a scalar meson (a) leads 
to a self-energy of the nucleon in the nuclear medium which contains a large scalar 
component S s and a large time like vector component S°. These two components 
compensate each other to a large extent if one calculates the single-particle energy. 
This leads to the well known fact that the binding-energy of nuclei are very small as 
compared to the rest mass of the nucleon M, a fact which has often been used to argue 
that relativistic effects should be small in the many-body problem of nuclear physics. 
The individual components of the self-energy S s are of the order of the nucleon rest 
mass (typically one third of M). Therefore the structure of the Dirac spinors in the 
nuclear medium is modified to quite some extent as compared to the Dirac spinor for a 
free nucleon. This medium dependence of the Dirac spinors, which affects the evalua- 
tion of the NN interaction in the medium, leads to a saturation mechanism for nuclear 
matter. 
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These more or less empirical models received support from calculations, which start 
from a realistic One-Boson- Exchange model of the NN interactional . In this context 
realistic NN interaction means that the parameters contained in these models are ad- 
justed to describe the experimental data of free NN scattering. Using such a realis- 
tic NN interaction in a many-body calculation of the Dirac-Brueckner-Hartree-Fock 
(DBHF) type, results were obtained for the saturation point of nuclear matter which 
were in quite a good agreement with the empirical saturation point [|3], M. Such a DBHF 
calculation accounts for the effect of correlations on the level of the BHF approxima- 
tion, i.e. the lowest order in the hole line expansion, and allows for the relativistic 
effects which we just outlined. This success of the relativistic features contained in the 
DBHF approach could be a solution of an old problem: the so-called Coester band 
phenomenon [||], which stands for the fact that many-body calculations based on vari- 
ous realistic models of the NN interaction lead to predictions for the saturation point 
of nuclear matter, which either fail to yield enough binding or predict a saturation 
density twice as large as the empirical value. 

Such a Coester band can also be observed in non-relativistic studies of finite nuc- 
lei @. Therefore it was quite obvious that attempts have been made to include the 
relativistic features also in DBHF calculations of finite nuclei. Indeed the inclusion of 
relativistic effects gave a substantial improvement in the calculated binding energies 
and radii of finite nuclei]?], |J. However, there remains still a discrepancy between the 
DBHF results and the experimental values. 

On the other hand it was observed that an extension of the non-relativistic BHF 
approach to a definition of the nucleon self-energy which accounts in a symmetric way 
not only for the particle-particle ladders included in the Brueckner G-matrix but also 
for the corresponding hole-hole scattering term may lead to a larger binding energy 
and a larger radius than obtained in the BHF approach ||. A similar feature, mov- 
ing the saturation point away from the Coester band by including hole-hole scattering 
terms, has also been observed in particle-particle hole-hole RPA calculation of nuclear 



matter [|10|, lUfl . Taking a very optimistic point of view one may argue that the com- 
bination of relativistic effects and hole-hole scattering terms may lead to an improved 
microscopic understanding of groundstate properties of nuclear systems. 

In the work presented here we would like to investigate the effects of particle-particle 
and hole-hole scattering terms on the relativistic structure of the nucleon self-energy. 
For that purpose we consider the nucleon self-energy defined with all terms of first and 
second order in the NN interaction (see Fig. 1). In contrast to earlier invest igat ions [12, 
TB| , |T3[ |T5| we evaluate the imaginary contributions to the nucleon self-energy allowing 



for all possible combinations of momentum k and energy A; . Therefore we can use 
dispersion relations to evaluate also the real part of the second order diagrams. We 
investigate in detail the effect of these higher order diagrams on the Dirac structure 
of the self-energy. Furthermore we derive from this relativistic self-energy an optical 
potential to be used in a Schrodinger equation for nucleon-nucleus scattering. 

The technical details for the calculation of the imaginary components in the self- 
energy are presented in section 2. The dispersion relations used to evaluate the corre- 
sponding real components are given in section 3. That section also contains a detailed 
discussion of the results. The main conclusions are summarized in the final section. 
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Figure 1: Graphical representation of the Hartree (a), Fock (b), 2plh direct (c), 2plh 
exchange (d), 2hlp direct (e) and 2hlp exchange (f) contributions to the self-energy of 
the nucleon. The momenta labelling the various contractions in (c) correspond to the 
nomenclature used in ([15] ) 



2 Imaginary part of the nucleon self-energy 

The physical system of neutron or nuclear matter is, by definition, translational as well 
as rotational invariant. Furthermore it is also assumed to be invariant under parity 
and time reversal. In general, under these conditions, the nucleon self-energy can be 
writtenP] as [|I|] : 

E(k) = E s (k) - 7o S°(k) + 7 • kE^(k) (1) 

where the functions E r (r = s, 0,v) can be projected out by taking the appropriate 
traces: 



E s 
E° 
E" 



^rr[ 7o E] 



(2) 
(3) 
(4) 



As mentioned in the introduction, in this work we want to go beyond the usual 
Hartree-Fock approximation to the self-energy, visualized by the diagrams a and b 
of Fig. 1, and pay special attention to the 2hole-lparticle (2hlp) and 2particle-lhole 
(2plh) diagrams displayed in Figs, lc, Id, le and If. Therefore, the different com- 
ponents of the self-energy are complex functions and, in general, we will write them 
as: 



V r + iW r 



(5) 



1 Notation: within formulae we shall use roman type (p, q, • • •) for 4- vectors, boldface type (p, q, • • •) 
for 3-vectors and normal math italic type (p, q, ■ ■ ■) for the norm of 3-vectors. 
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Previous works have focused on the imaginary part of the self-energy at the on-shell 
energy of the propagating nucleon. In Ref. [0] the 2plh and 2hlp direct contributions 



where studied, while the work of Ref. fbl considered also the Fock-exchange terms. 



This study, however, was restricted to the on-shell 2hlp contributions. 

In the present work we extend these calculations and evaluate the direct and ex- 
change terms for both 2hlp and 2plh contributions to the nucleon self-energy consid- 
ering off-shell effects, i.e. investigating the self-energy E(&o,k) for all combinations of 
energy ko and momentum k. In a first step we calculate the imaginary part of the 
self-energy, as described below, from which we later obtain the real part by means of 
a dispersion relation. 

The nucleon-nucleon interaction is derived from a G-matrix evaluated within the 
Dirac-Brueckner-Hartree-Fock approach. This G-matrix is parameterized in terms of 
an exchange of effective a and u mesons|7|, 16|. The starting point of our discussion 



will be the Hartree-Fock (HF) approach. Therefore, and for the sake of simplifying the 
notation, we will define the following effective quantities in terms of the Hartree-Fock 
components of the self-energy: 

M*(k) = M + E s HF (k) (6) 
k* = k (l + S% F (k )) (7) 
E*(k) = Vk* 2 + M* 2 (8) 

k* = (k*,k*) = (k + X° HF (k),k*) (9) 
e HF (k) = E*(k)-X° HF (k) (10) 

The nucleon HF propagator is then 
where 

= 6(E*(k)-E F ) d(E F -E*(k)) _ 1 {m 

S{ ' kl-E*{k)+iri k%-E*(k)-ir] k% + E*(k)-ir} { ' 

with E* F = <J k* F + M* 2 . In eq. ( |12"D we will disregard the last term, which is the 
anti-nucleon contribution. This is quite reasonable due to two reasons. First of all, the 
production of fermion anti-fermion pairs should be negligible at the energy scales of 
interest and secondly the pair production is inhibited in high density matter ]T2"| , |TB|] . 
In the actual calculations we set Yj^jjp — and disregard the momentum dependence 
of an d ^°hf f° r which we take their values at the Fermi momentum kp. This 
approximation relies on the weak momentum dependence of the Hartree-Fock self- 
energies [0. In this way, M* acquires a fixed value, k* = k, and the only effective 
quantity that depends on the nucleon momentum is E^, which from now on will be 
labelled with a subindex. It is also useful to split up the fermion-propagator into 
a particle and a hole term and consider explicitly all terms of Fig. 1 by taking the 
following prescriptions for the nucleon propagators in the diagrams: 

V* _|_ ]\/J* 

Particle line : iG p (k) = i — — — g p (k) (13) 

Hole line : iG h (k) = i ^- g h (k) , (14) 

2E k 
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where the (^-functions for particle and hole states are defined as: 

0{K - E F ) 



9h<$ 



k* -E£ + i V 
9(E* F - Ei) 



These ^-functions are formally similar to the non-relativistic HF nucleon propagators. 
Therefore, applying the standard Feynman rules together with the prescriptions of Eqs. 
(|13|) and (14), we obtain a general expression for the 2plh contribution 

r d 4 q d 4 p* 

xL x . ab (kQ, k) = ij — — A M X ;ab% P, q)^ P (P + q)^(p)^ P (k - q) (15) 

while the formula for the 2hlp state is identical to the above except for the exchange 
of the p and h subscripts in the ^-functions. The subscript X stands for D (direct) 
and E (exchange) diagrams, while a and b stand for the different meson types. 

The function M. x - a b{k, p, q) contains, basically, the Dirac structure of the interac- 
tion and, for direct terms, it has the general form, 

\ 2 2 

M D;ah (k,p,q) = ~' 9a % r 6 (q) A 6 (q) rr[r 6 (q) (f+ (f + M*) r a (q) 



Q 77"* 771* 77 1 * 



(f + M*)](Jf- ff + M*) A a (q) r a (q) (16) 



The constant Xj is the isospin degeneracy and comes from the loop trace, g a and g b are 
the meson-nucleon coupling constants. The meson-nucleon vertices are denoted by V 
and meson propagators by A. A form factor of the type F(q) = A £i_ 2 , with a typical 
cut-off mass of A = 1500 MeV, has been attached to each vertex. This is equivalent to 
modifying the meson propagators in the following way: 

A a (q)^A a (q)F 2 (q) (17) 

Note as well the minus sign that comes from the fermion loop. For exchange terms one 
can write: 

M E]ab (K P, q) = ^S 9 K* T ^ ~ <* ~ ^ - * ~ P) ^ + M *) 

8E p E p+q E k _ g 

r.(q) (f + M*) r 6 (k - q - p) (K*- Jf + M*) A a (q) T a (q) (18) 
The imaginary part of the self-energy is obtained as: 

W x ., ab (k ,k) = 4vr 3 j ^^M x , ab (k,p,q) 6(k,p,q) (19) 

Where, for 2plh, the factor has the following form: 

6(k, p, q) = 6(Ep - p*)9(k* -q - E F )6(p* + q - E* F ) 

Kvt - E;)5(p* + q - E; +q )5(k* -q - E*_ q ) (20) 



and, for 2hlp: 

0(k, p, q) = -6(p* - E* F )9(E F -k* + q )9(E F -p* - q ) 

5(p* - E;)6( P * q + g - E; +q )5(k* - q - £*_ q ) (21) 

One can now use the step functions to define the boundaries for the integrations over 
dqo and dp - Working in spherical coordinates we can automatically perform one of 
the axial integrations (say over the axial q- angle). By rewriting the first 5-function in 
terms of the 3-momentum variable, the dp integration can be easily performed. After 
that we obtain, for the 2plh diagrams, 



E% 



W x -ab(ko,k) = 4 J dq J dq q 2 d(cos6 q ) J dp* p J dip p d(cos9 p ) 



2(2tt) / 

M X;ab (k, p, q) E* p S(p* + q - E; +q )5(k* - g - E*_ q ) (22) 



PO 



Where, p takes the value p = yPo 2 — M* 2 , and 

Po = max{E F -qo,M*). (23) 

For the 2hlp contribution one obtains a similar expression up to an overall minus sign. 
The polar integration over the angle 9 q between q and the external momentum k (taken 
along the z-axis) is readily performed by using the delta function 



r/ t * F * ,_K- qx( a 2A:ggo-q 2 + M* 2 -k . 
6{k - g - £ k _ q ) = 5(cos 6 kq — ) , 

which, in order to keep | cos 6k q | < 1, imposes the following constraints to the q variable: 

q- <q<q+ 

where: 

g ± = b±v/(A:o*-go) 2 -M* 2 | . 

The integration over the polar angle 9 p can also be performed very easily by referring 
the two angular variables of the momentum p (4>p,9 p ) to a reference frame in which q 
acts as the z-axis and using the remaining 5-function 

5(p Q + g - E p+q ) = 5(cos 9 pq — ) . 

The requirement that the absolute value of the cosine must be less than one puts 
restrictions to the values of Pq, similar to those obtained for g, but more complicated 
to be implemented analytically due to the two different sign possibilities for q 2 . For this 
reason, we have taken care of the restrictions over p^ numerically through an explicit 
step function inside the integral. 

After all these considerations we can write our final expression for the imaginary 
part of the self-energy. For the 2plh state we have: 



r *2 



K~ E F 1+ E F 



2(2n) 4 k 

q- po 



-u —t- ^+ - 1- [n * I 2"| 2 2lT 



2pq 

J o 

Mx;ab(K' k iPoiqO,q,<Pp), (24) 
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Figure 2: Imaginary part of the self-energy calculated for a fixed momentum k = kp 
= 1.4 fm _1 as a function of the energy variable, which is normalized such that it is 
zero at the Fermi energy. Results are displayed for the scalar part (W s , solid line), 
the time-like vector part (W°, dashed line) and the space-like vector component (W v , 
dashed dotted line) 



and for the 2hlp state: 

W X ;ab{k , k) = . , 4 J dq J dq J dpi 6(1 - 

k*-E* F q- Ep 

Mx;ab(koi k iPoi<lo,q,iPp) , (25) 

where it is understood that M. does not contain the energies in the denominators 
of Eqs. ( |i~6|) and (|18j) becasue they have canceled out with the energy factors that 
appeared on rewriting the 5-functions. 



2pggo + q : 
2pq 



■ 1 * 



dip F 



3 Results and Discussion 

In the first part of this section we will concentrate the discussion of results on the 
example of nuclear matter at a density close to the saturation density, in particular we 
will consider a Fermi momentum kp of 1.4 fm _1 . We are going to discuss the contribu- 
tion of the diagrams of second order in the NN interaction to the various components 
of the imaginary part in the nucleon self-energy. The NN interaction is described in 
terms of the exchange of an effective scalar, a, meson and an effective vector, to, me- 
son. The masses of these mesons are fixed to the masses of the corresponding mesons 
in a realistic meson exchange model of the NN interaction (m a = 550 MeV, m u = 
782.6 MeV), while the effective coupling constants have been adjusted such that 
a Dirac-Hartree-Fock calculation of nuclear matter at kp = 1.4 fm _1 using these con- 
stants would reproduce the results for single-particle energies and the binding energy 
obtained in a Dirac-Brueckner-Hartree-Fock (DBHF) calculation employing Bonn A 
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NN potential^]. This Dirac-Hartree-Fock parameterization of the DBHF results yields 
at this density coupling constants of ^=8.536 and ^=9.536 for the a and the u) meson, 
respectively 0. The single-particle Green function G of eq. (|TT|) for the propagation of 
the intermediate states is also defined employing the results of this DBHF calculation. 
This means that we consider a HF self-energy defined in terms of 

E S HF = -374.9 MeV 

E° HF = -289.8 MeV. (26) 

As a first example we show in Fig. 2 the imaginary part of the self-energy calculated 
for nucleons with a momentum fixed to k = kp = 1.4 fm _1 as a function of the energy. 
Note that the energy variable in this figure is normalized such that an energy zero 
corresponds to the Fermi energy ep = ^HF(kp) as defined in (|ItJ). Results are displayed 
for the scalar component, ImagS s = W s , the time-like vector component, ImagS = 
W°, and the space-like vector component, klmagE v = W v . All these components are 
negative at energies below the Fermi energy and positive for those above. They are 
much larger at positive energies, reflecting the fact that the phase space of 2 particle 1 
hole (2plh) states with this momentum k is considerably larger than the corresponding 
phase space of 2 hole 1 particle states (2hlp). The scalar and time-like vector part are 
of similar size and exhibit the same sign. Using our notation this means that these 
contributions cancel each other to a large extent in calculating the expectation value 
for a Dirac spinor u representing a nucleon, i.e. a solution of the Dirac equation at 
positive energy: 

M* kk* 
«(k)ImagS(A;o, k)it(k) = -W°(k , k) + -^W s (k , k) + — (27) 

using the definition 

k = kW v (k ,k) (28) 

and E*, M*, k* as defined in eqs.@ - (|8]). The Dirac spinors u are normalized such 
that u'u = 1. As the absolute value of W° is always larger than the absolute value 
of W s at the same energy, this expectation value, which is roughly proportional to 
W s — W°, shall be positive at energies below the Fermi energy and negative above, 
as it is the case for the imaginary part of the self-energy calculated within a non- 
relativistic frame. The space-component of the vector part, k, is significantly smaller 
than the other terms. This difference, however, is not as large as one finds, e.g. in the 
real part of the self-energy calculated within the Dirac-Hartree-Fock approach. Similar 
results are obtained for other momenta and nuclear densities. 

From the 2hlp contribution to the imaginary part, i.e. the one at energies below 
the Fermi energy E* F one can determine the corresponding contribution to the real part 
by applying the dispersion relation 

AV& lp (u,k) = RealAE^u, k) 

= Pro ^hnag^Kj) 

where the P is used to indicate the principle value integral and the index a represents 
s, 0, and v, referring to the scalar and vector components of S. Note that here and in 
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Figure 3: Real part of the 2hlp contribution to the self-energy for a fixed momentum 
as a function of energy. For further details see Fig. 2. 

the following the energy variables u are redefined such that u = corresponds to the 
Fermi energy k = ep. Typical examples for the real part of the 2hlp contribution to the 
self-energy are displayed in Fig. 3. The energy dependence of these terms, can easily be 
understood from the energy dependence of the imaginary part, shown in Fig. 2, and the 
dispersion relation of eq. (|29"D . All three contributions are negative for energies between 
—100 MeV and the Fermi energy, which are typical energies for quasihole states, as 
well as for positive energies, i.e. particle states with energies above the Fermi energy 
e F . As the absolute value of the time-like vector component RealAE^ (dashed line) 
is consistently larger than the corresponding scalar component (solid line), we obtain 
a repulsive contribution to the quasiparticle energy, arising from this 2hlp term from 
energies starting around —100 MeV below the Fermi energy up to infinity . This 
repulsive contribution is largest for quasihole states with energies below €p and will 
decrease for energies above ep with increasing energy. It should be noted that the 
space-like vector component exhibits a sizable contribution. 

The real part of the 2plh contribution to the nucleon self-energy can be calculated 
from the corresponding imaginary part by a dispersion relation rather similar to eq.(|2"§D 



P r°° ,ImagESL 1 i,(a/, k) 
Rea\AE% lh (u,k) = --J o duo' ; . (30) 

We already discussed above that the 2plh states lead to larger contributions to the 
imaginary part than the 2hlp terms. Therefore it is clear that the 2plh contributions 
to the real part of the self-energy are significantly larger than those originating from 
2hlp terms. This can be seen from Fig. 4, which employs a scale which is about a 
factor 10 larger than the one used in Fig. 3 to visualize the 2hlp contributions. 

In the energy range of interest all three terms are positive. This means that the 2plh 
contribution to the scalar part of the self-energy E s tends to compensate the negative 
Hartree-Fock contribution to this term originating mainly from the exchange of two 
correlated pions, which is parameterized in the realistic OBE potentials by means of 
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Figure 4: Real part of the 2plh contribution to the self-energy for a fixed momentum 
as a function of energy. For further details see Fig. 2. 

the a meson exchange. A similar situation also arises in the case of the timelike vector 
component of the nucleon self-energy S°: The negative Hartree-Fock contribution to 
E°, which is mainly due to the exchange of the to meson, is compensated to some extent 
by the 2plh terms. This means that the 2plh terms tend to reduce the Hartree-Fock 
contributions to the various terms in the self-energy while the 2hlp corrections yield 
contributions to S s and S° with the same sign as the Hartree-Fock terms. Similar 
results are obtained for other momenta k. 

Since, however, our Hartree-Fock approximation to the self-energy has already been 
extracted from a DBHF calculation, we are not allowed to simply add the real part of 
2plh contribution to the self-energy to the corresponding DBHF results. This would 
lead to a double counting of these 2plh terms. Instead we use a subtracted dispersion 
relation defined by 

&V% lh (u, k) = RealA£« k) - RealAE% lh (e HF (k) - e F , k) (31) 

using the real contributions to the various terms a in the self-energy as defined in 
(p0|) . With these definitions we now define a quasiparticle self-energy, which is real and 
energy dependent, to be 

V«{u, k) = Z a HF (k) + AV 2 a plh (u, k) + AV^>, k) (32) 

and we can determine a quasiparticle spinor 



u qp (k) 

with 



E* + M* 1 

qp ' qp I 



\ 2E* 



M* qp = M + V* p (u qp ,k) 
k; p = k(l + V q v p (u qp ,k) 
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Figure 5: Results for the quasiparticle energy, calculated according to eq. fl3"E|), are 
compared to the correponding energies calculated in the HF approximation for the 
self-energy (dashed line) 



E* qp = ^) 2 + (m;^ 

o 



UJ, 



E* ~ V q Uu qp , k)-e F . (34) 



qp qp qp 



This means that this spinor is an eigenstate of the Dirac equation 

1 + V q v p )a ■ k + 7 °(M + V q %) - V q ° p ] u qp = e qp u qp , (35) 



which uses the self-energy V qp of (|32| ) calculated at the energy uj which corresponds to 
the quasiparticle energy defined by 

U*0 = E* qp - El(u qp , k) . (36) 

If for the moment we ignore the 2hlp contribution, AV^ lp , to the quasiparticle self- 
energy in (|32|) , the subtraction defined in (31) ensures that the AV^ lh terms vanish on- 
shell and the quasiparticle energy coincides with the HF energy, thus avoiding double 
counting. However, when the contribution of the 2hlp terms is taken into account, 
the self-consistent definition of the energy variable u qp gives rise to a non-vanishing 
correction due to energy dependence of the 2plh terms. 

Results for the quasiparticle energy are displayed in Fig. 5 for various momenta k. 
It can bee seen that the inclusion of the 2hlp terms yields a significant reduction in the 
quasiparticle energy e qp as compared to the corresponding HF result. This reduction 
is as large as 17.5 MeV for momenta close to zero and reduces to values around 5.3 
MeV for k = kp and becomes negligibly small around k = Ikp. This repulsive effect 
of the 2hlp terms in the quasiparticle energy has been discussed already in context 
with Fig. 3. Such a large effect for the quasiparticle energy may indicate that the 
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2hlp terms may as well have some effect in the calculation of the total energy. For the 
calculation of the total energy, however, it is not sufficient to evaluate the quasiparticle 
energy, but one would need the whole spectral distribution of hole strength for states 



with momenta below and above the Fermi momentum |T7[ . 

It is one aim of our study to explore the effects of the Dirac structure of the nucleon 
self-energy, calculated beyond the mean field approximation, on the complex optical 
potential for nucleon-nucleus scattering. The Dirac equation reads now: 

1 + E^)« • k + 7 °(M + E* p ) - EjJ u qp = ~e qp u qp (37) 



where the full complex self-energy is used. Real e qp solutions of eq. (|37j) involve complex 
values of k. We take as approximate solutions of eq. fl3~7|) the values e qp determined 
from eq. (|36l). It is convenient to rewrite eq. ( |37|) into a form which only contains an 
effective scalar V s and vector potential V [P§ 



k + 7 °(M + V s ) - V° 



U qp tqpUqp 



where 

V s 



V° = qp 9P qp (38) 
1 + E«, 

The self-energy terms T, qp have been calculated in the quasiparticle approach of (|32| ) us- 
ing the self-consistent relation between the three-momentum and the energy as defined 
in ([TJ]). This Dirac equation can be transformed into a Schrodinger- type equation for 
the large component of the Dirac spinor, leading to a complex and energy- dependent 
Schrodinger equivalent potential of the form 

u = ys ~ + m [ (ys)2 " (V ° )2 1 • (39) 

Results for the real part of the renormalized components V s and V° are displayed 
in Fig. 6 as a function of the energy variable u qp , i.e. the quasiparticle energy e qp 
normalized such that the Fermi energy occurs at zero. For a comparison we also show 
the unrenormalized components V q a p . The difference between the solid and dashed line 
is a measure for the importance of the space-like vector component E u of the self-energy. 
We find that this space like components yield a slight reduction of the absolute values 
for RealV, which is of the order of 3 percent. Larger effects only occur at negative 
energies, below the Fermi energy, where the 2hlp contribution gets more important. 

This figure also demonstrates that the energy- dependence of the various Dirac com- 
ponents in the real part of the self-energy remains weak, again with the exception of 
energies below ep. Also the deviations from the Hartree-Fock result, -375 MeV and 
-290 MeV for the scalar and vector parts, respectively, are not very pronounced. This 
can also be seen from the left part of Fig. 7, exhibiting the real part of the Schrodinger 
equivalent potential U defined in ([53). The results derived from the quasiparticle ap- 
proximation (solid line) including the effects of the energy-dependence in the 2plh and 
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Figure 6: The real parts of the bare scalar and vector components (T, s qp and S° p ), rep- 
resented by the dashed line, are compared to the real parts of V s and V°, 
according to eq. (pq). 



renormalized 




Figure 7: Energy dependence of the real (left) and imaginary part (right) of the 
Schrodinger equivalent potential. The result for the real part is compared to the 
corresponding prediction obtained in the HF approximation. The dashed curve in the 
right part of the figure exhibits the expectation value of W calculated according to 

<m>. 
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Figure 8: Bare and renormalized scalar and vector components of the imaginary part 
of the quasiparticle self-energy. For further details see Fig. 6. 



2hlp terms, exhibit a dependence on the quasiparticle energy, which is very similar 
to the one derived from the HF approximation. Therefore one may conclude that the 
energy-dependence in the depth of central, Woods-Saxon type, optical potentials, used 
to describe nucleon-nucleus scattering is mainly due to the relativistic structure of the 
underlying Hartree-Fock self-energy to be used in a Dirac equation. Dispersive effects 
due to the consideration of 2plh and 2hlp terms also lead to an energy dependence, 
which, however, is much smaller. These dispersive corrections tend to make the po- 
tential slightly more attractive at higher energies. It should be mentioned, however, 
that an energy dependence of the central Schrodinger potential similar to the one ob- 
tained here within the relativistic scheme can also be obtained within a non-relativistic 



Hartree-Fock due to large non-localities in the Hartree-Fock potential HT9 



An analysis rather similar to the one just outlined for the real components of the 
self-energy can also be performed for the imaginary parts. Results for the renormalized 
Dirac components W a = ImagV" (a = s and 0) are presented in Fig. 8 and compared 
to the bare scalar- and vector imaginary components. As we discussed already in 
the beginning of this section, the imaginary part of the bare components are negative 
for energies below the Fermi energy (2hlp), zero around u = (ep) and positive for 
positive energies (2plh). The renormalizing effects of the space like vector components 
are much more important for the imaginary part than for the real part. This can be 
deduced from the differences between W a and W a : the bare components are as large 
as twice the renormalized quantities. 

The imaginary part of the Schrodinger optical potential [eq. (p9[)j is shown on 
the right of Fig. 7. An alternative way of evaluating this imaginary part would be to 
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Figure 9: Real and Imaginary parts of the Schrodinger equivalent potential, calculated 
using the full model (solid lines), ignoring exchange diagrams in the evaluation of the 
terms of second order (dashed line), and using the Hartree approximation also for the 
lowest order term (dashed-dotted line) 



calculate the expectation value of the Dirac operator 

W s - j°W° + 7 • kW v ] u qp (40) 



using the quasiparticle Dirac spinors of fl33|). These expectation values lead to the 
dashed line on the right part of Fig. 7. We see that the results for the imaginary part 
of the Schrodinger potential are rather insensitive on the way of calculation. 

In the calculation of the imaginary components of the self-energy discussed in sec- 
tion 2 and consequently also in the corresponding real components, we always deter- 
mined the Direct- and Exchange-contributions (see eqs. (^) and (0)). In order to 
explore the importance of the 2plh and 2hlp exchange terms we compare in Fig. 9 
the results obtained for the real and imaginary part of the Schrodinger equivalent po- 
tential with (solid lines) and without (dashed lines) inclusion of the exchange terms in 
the 2plh and 2hlp parts of the self-energy. One finds that the effects of the exchange 
terms on the real part are rather small typically around 1 MeV. However, the effects 
are significantly larger for the imaginary part where the difference gets as large as 30 
percent of the total result. This can be understood from the fact that the real part 
of the self-energy is dominated by the Hartree-Fock contribution, which is identical in 
these two approaches. 

If one ignores the exchange terms in calculating the 2plh and 2hlp terms, one may 
consider it more consistent to ignore the exchange term also in the leading contribution 
and replace the Hartree-Fock approximation by the Hartree-approach. In ref. J7| effec- 
tive meson-nucleon coupling constants were determined to reproduce the DBHF results 
within a Dirac-Hartree model. The resulting coupling constants are a bit larger than 
those derived from the Dirac-Hartree-Fock analysis. If we use these Hartree-coupling 
constants and ignore the effects of exchange terms in the leading term as well as in 
the 2plh and 2hlp terms, one arrives at a Schrodinger equivalent potential displayed 
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Figure 10: Real and Imaginary parts of the Schrodinger equivalent potential calculated 
for nuclear matter with a Fermi momentum of kp = 1.4 fm -1 (solid line) are compared 
to those obtained for &p = 1.2 fm -1 . 

by the dashed-dotted lines in Fig. 9. The differences to the solid lines get now quite 
pronounced for the real part. This can be traced back to the fact that in the Dirac- 
Hartree approximation S s and S° are independent on the momentum or energy of the 
state and T, v is identical to zero. The results for the imaginary part obtained in this 
approach, however, are rather close to those evaluated with inclusion of the echange 
contributions. 

It is also interesting to investigate the density dependence of the self-energy in order 
to extend these calculations to finite nuclei. As an example we present some results 
obtained for nuclear matter with a Fermi momentum kp = 1.2 fm -1 in Fig. 10. The 
results obtained at these various densities, either for the self-energy keeping track of 
the Dirac structure, or using the Schrodinger equivalent potentials derived from these 
components, may than be used in a local density approximation for a prediction of 
nucleon-nucleus scattering. Such investigations are in progress. 

4 Conclusions 

The relativistic structure of the self-energy for a nucleon in nuclear matter is investi- 
gated by including all irreducible terms of first and second order in the residual interac- 
tion. For the NN interaction a parameterization of the Dirac-Brueckner-Hartree-Fock 
(DBHF) G-matrix in terms of the exchange of effective scalar and vector mesons has 
been used. The 2plh and 2hlp contributions to the imaginary part of the self-energy 
ar e evaluated keeping track of all direct and exchange terms. The corresponding 2plh 
and 2hlp contributions to the real part are derived from these imaginary components 
by means of dispersion relations. A subtracted dispersion relation must be used for the 
2plh term to avoid double-counting with the G- matrix underlying the DBHF approach 
on which these studies are based. 

The inclusion of 2hlp diagrams in the evaluation of the real part of the self-energy 
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yields a non-negligible modification of the scalar and vector components in particular 
for states with momenta below the Fermi momentum. Also the value of the quasipar- 
ticle energy is increased by a value as large as 17.5 MeV for momenta close to zero 
and to values around 5.3 MeV for k = k F . Such a large effect for the quasiparticle 
energy may indicate that the 2hlp terms should have some effect in the calculation of 
the total energy. 

The calculated self-energy can also be transformed into a Schrodinger equivalent 
optical potential, to be used in the study of nucleon-nucleus scattering. Exchange 
diagrams are no n- negligible in the evaluation of the imaginary part. The energy- or 
momentum-dependence of the central component of the real potential, however, is 
dominated by the effects of the Dirac-Hartree-Fock contribution. The 2plh and 2hlp 
terms give rise to a more attractive SEP at positive enrgies and introduce an additional 
energy dependence which is very weak. 
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